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We demonstrate that the recent measurement of azimuthally collimated, long-range rapidity 
(“ridge”) correlations in ^/s = 13 TeV proton-proton (p-fp) collisions by the ATLAS collabora¬ 
tion at the LHC are in agreement with expectations from the color glass condensate effective theory 
of high energy QCD. The observation that the integrated near side yield as a function of multiplicity 
is independent of collision energy is a natural consequence of the fact that multiparticle production 
is driven by a single semihard saturation scale in the color glass condensate framework. We argue 
further that the azimuthal structure of these recent ATLAS ridge measurements strongly constrain 
hydrodynamic interpretations of such correlations in high multiplicity p+p collisions. 


The systematics of azimuthally correlated, long-range 
rapidity correlations-the so-called ridge correlations- 
have recently been measured in small systems rang¬ 
ing from proton-proton (p+p), proton-nucleus (p+A), 
deuteron-nucleus (d+A) and Helium 3-nucleus (^He+A) 
collisions at both RHIC and LHC HHT]. An outstand¬ 
ing question is whether such correlations arise from the 
final state hydrodynamic flow of the produced matter- 
as is believed to be the case for nucleus-nucleus (A+A) 
collisions-or from correlations that are already embed¬ 
ded in the initial state wavefunctions of the projectile 
and target. 

In a series of papers, two of us argued that many 
features of the ridge correlations in high multiplicity 
p+p collisions could be quantitatively explained within 
the initial state framework of the color glass condensate 
(CGC) effective theory of high energy QCD [SUIT]. Our 
results were derived within a “glasma graph” approxi¬ 
mation to the most general set of possible two particle 
correlation Feynman diagrams that contribute within the 
power counting scheme established by the CGC effective 
theory. 

Subsequent to the early p+p high multiplicity ridge 
measurements, ridge correlations observed in light-heavy 
ion collisions at both the RHIC and LHC appear to indi¬ 
cate a collective origin of some of the systematics exhib¬ 
ited by the measurements [T2UT^ . However some of these 
systematics adduced to support the hydrodynamic flow 
interpretation may also arise within initial state frame¬ 
works [TM2H] and others, notably the hadron fragmen¬ 
tation patterns often attributed to hydrodynamical flow, 
may have alternative interpretations PSI150] . We shall 
not discuss these further here but merely note that the 
interpretations of the origins of the ridge correlations re¬ 
main fluid; further data is urgently needed to help sort 
through the competing (and compelling) explanations. 

Recently the ATLAS Collaboration reported first mea¬ 
surements of ridge correlations in high multiplicity p+p 
collisions at the LHC with the highest center-of-mass en¬ 
ergy of 13 TeV. They found quite strikingly that the as¬ 
sociated integrated yield of the ridge collimation agrees 
(within experimental errors) with the CMS results at the 


lower 7 TeV center-of-mass energy when plotted as a 
function of multiplicity |31j . We will demonstrate here 
that the observed energy independence of the ridge is a 
natural consequence of the fact that multiparticle pro¬ 
duction in the color glass condensate is controlled by a 
single semihard saturation scale Qs in collisions of sym¬ 
metric systems like p+p. More concretely, we will show 
that the ATLAS results are in quantitative agreement 
with our prior computations with no adjustment of pa¬ 
rameters [55] . 

Before we proceed to discuss the details of the glasma 
graph computation, it is useful to keep in mind a few 
generic features of multiparticle production in this frame¬ 
work. The azimuthal structure of the two particle corre¬ 
lations results from an interplay between mini-jet (more 
generically di-jet) production along with glasma graph 
contributions at high gluon densities ISIIIQ]. 

Di-jet production proceeds via single gluon exchange in 
the t-channel and is the dominant process for di-hadron 
production in minimum bias events corresponding to low 
parton densities. As the two gluons are produced from 
the same t-channel gluon exchange, the correlation is 
kinematically constrained to be back-to-back in their rel¬ 
ative azimuthal angle (A$ « tt). In contrast, the glasma 
graphs contain two independent ladders and are there¬ 
fore not similarly constrained; they are however highly 
suppressed for high transverse momenta and low parton 
densities. It was argued however [551|51j, based on the 
power counting of the GGC effective theory, that glasma 
graphs can be enhanced by with respect to the di¬ 
jet contributions at high parton densities. The fact that 
the proton unintegrated gluon distributions are peaked 
about transverse momenta Qs in the CGC effective the¬ 
ory, and the kinematical constraints that the two gluons 
are produced with momenta pp ^ <lT ^ Qs generates an¬ 
gular correlations that are symmetric around A$ = -k 12 
and enhanced at A$ « 0 and A$ « tt. The glasma 
graph two particle correlations therefore provide a natu¬ 
ral source of the symmetric near side and away side ridge 
in high multiplicity collisions. 

We note further that, because of the azimuthal symme¬ 
try noted, only positive even Fourier coefficients |35j of 
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the azimuthal dependence of the two particle glasma cor¬ 
relation can be entertained in the glasma graph approx¬ 
imation. On the other hand, the Fourier decomposition 
of the di-jet gives both even and odd Fourier coefficients. 
The restriction of the glasma correlation to even Fourier 
coefficients is a consequence of the glasma graph approx¬ 
imation which is valid for but Qs ^ PttQt^ 

where Aqcd 100 — 200 MeV is the intrinsic QCD scale. 
For pt^Qt ^ Qs, positive odd moments in the Fourier de¬ 
composition of the two gluon azimuthal correlation are 
also feasible [55J [35] due to multiple scattering contribu¬ 
tions that break the charge and parity invariance mani¬ 
fest in the glasma graph computation. Thus depending 
on the values of Qs and pt,<It examined, an interest¬ 
ing pattern of positive and negative Fourier modes of the 
azimuthal correlations may be anticipated. In contrast, 
hydrodynamics will generically give Fourier coefhcients of 
the two particle correlation that are positive definite m- 
We will return to this point later on in the manuscript. 

Last but not least, the glasma framework is straight¬ 
forwardly extended to three particle correlations [38| and 
generalized to n-particle correlations as well [33]. In the 
latter case, the resulting cumulants naturally result in 
a multiplicity distribution having the form of a negative 
binomial distribution. The tail of this distribution cor¬ 
responds to the rare high multiplicity events relevant for 
our discussion This point too shall be revisited be¬ 
fore we conclude our discussion. 

Now turning to the glasma graph computation, we first 
note that the double inclusive gluon distribution is ap¬ 
proximated by a fer-factorized form m- There are a 
total of eight glasma graphs that are used in this analy¬ 
sis; the full expressions can be found in m- Let us give 
as an example the expression from one such diagram re¬ 
sponsible for the near side ridge 
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where the unintegrated gluon distributions (UGDs) <i> 
are evaluated at (small) x values on the order x ^ 
Pt/- the precise prescription is given in the ref¬ 
erence above. The unintegrated gluon distributions can 
be expressed in terms of T{x, rj_)-the forward scattering 
amplitude of a quark-antiquark dipole of transverse size 
r_L on a proton/nuclear target-through the expression 
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The unintegrated gluon distribution at a given momen¬ 
tum fraction x is obtained by solving the Balitsky- 
Kovchegov (BK) renormalization equation |42|[43]. The 
implementation we use for our study includes the run¬ 
ning coupling next-to-leading-log (NLL) correction to 


BK, hereafter referred as the rcBK equation [33]. The 
initial condition for TA,Bix,r±) in the rcBK equation 
is set according to the McLerran-Venugopalan (MV) 
model |451 146] with a finite anomalous dimension. In 
the MV model, the shape of TA,B{x,r±) is expressed as 
a function of the dimensionless quantity r^Qo, where Qq 
is a nonperturbative scale. Fixing the value of Qq pre¬ 
scribes the initial shape of Ta,b{x, r_L) at an initial x = xo. 
However such an initialization is not unique since the 
scattering amplitude is known to have a strong depen¬ 
dence on the impact parameter m- Such a dependence 
can be accounted for by varying the value of Q§ to ini¬ 
tialize 7A,B{x,r±). 

The value of Qq — 0.168 GeV^ (in the fundamen¬ 
tal representation) describes the minimum bias (37] 135] 
HERA deeply inelastic electron-proton scattering data. 
Therefore higher Qq will correspond to the scattering of 
a color dipole off the proton at smaller impact param¬ 
eters where the color charge densities in the target are 
higher. Even if a fixed impact parameter is probed by 
the dipole, there could be fluctuations in the configura¬ 
tions of the color charge inside the proton purely driven 
by stochastic processes that are nonperturbative at the 
initial scale [32 ED] . 

The effect of such fluctuations will give rise to a distri¬ 
bution of Q)’s. By choosing higher Q§’s, one can there¬ 
fore incorporate the physics of rare configurations of color 
charge densities mm- We will therefore model rare 
high multiplicity p-l-p collisions by using UGDs corre¬ 
sponding to higher values of Qq. It is important to note 
here that the expressions Eqj^ and Eqj^ are inclusive 
quantities and, strictly speaking, high multiplicity events 
do not necessarily correspond to using a larger initial Qo 
in the inclusive yield. In a fully first principles approach, 
in order to get the multiplicity distribution right, one 
must include impact parameter fluctuations, color charge 
fluctuations (which generate a negative-binomial distri¬ 
bution) and intrinsic gluon fluctuations of the initial sat¬ 
uration momentum (See |5l] for more details). Such a 
first principles treatment, combined with small-x evolu¬ 
tion (here included via running coupling BK), is currently 
out of reach in high energy QGD. We therefore resort to 
adjusting the initial saturation scale Qq to the match the 
multiplicity of interest. 

However once this initial condition is set, the quan¬ 
tum evolution in x is uniquely given by the rcBK equa¬ 
tion which evolves Ta,b{x, rj_) (and therefore the UGDs) 
to smaller values of x. The peak of the distribution 
^A.B{x,k±_) corresponds to the saturation scale Q^ at 
that X. For the typical kinematics probed by these exper¬ 
iments, we find saturation scales (in the adjoint represen¬ 
tation) on the order of Q^ ~ 1.5 (3.0) GeV for min-bias 
(central) collisions. It must be noted that the UGDs that 
go into the calculation of the two particle correlation do 
not include any free parameters and are well constrained 
by HERA data. 

In addition to the just discussed UGDs, Eq.[2 depends 
on the transverse overlap area S± and the nonperturba- 
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tive constant C- The latter estimates contributions to 
multiparticle production below the scale Qj. The trans¬ 
verse overlap area S± cancels out in the expression for 
the per-trigger yield. This leaves only the parameter C, 
which was extracted previously by an independent 
factorization analysis of data on n-particle multiplicity 
distributions |S11|33] and estimated to be C = 1/6- Nu¬ 
merical simulations by solving classical Yang-Mills also 
constrain the value of C to be of similar order [M] . 

We should emphasize that precisely this value was used 
in our prior fits to the p-bp CMS 7 TeV data. 

The di-jet contribution in this framework is estimated 

[551155] to be 
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where Cp = — l)/2Nc and Q is the BFKL Green’s 

function that generates gluon emissions between the glu¬ 
ons that fragment into triggered hadrons. The detailed 
NLL expression of Q can be found in Ref m- As previ¬ 
ously, the transverse overlap area S± in Eq. [^ cancels in 
the final expression of per-trigger di-hadron yield. 

The combined effect of the glasma graph and di-jet 
contributions to the total two particle correlation is given 
by 
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The decomposition of the above double-inclusive distri¬ 
bution into a sum of glasma and BFKL (di-jet) contri¬ 
butions is far from trivial. The di-jet diagram (BFKL) is 
a one-loop correction (and therefore as suppressed but 
enhanced by A)?) with respect to the Glasma diagrams. 
Its full calculation in the “dense-dense” limit of high mul¬ 
tiplicity collisions requires knowledge of the background- 
field gluon propagator in these collisions, which is not yet 
available. 

Nevertheless, significant work has been done in the 
dilute dense limit of the CGG m 158] . where one ob¬ 
tains a contribution from the Glasma graph, the mini-jet 
(BFKL) graph as well as a contribution from the inter¬ 
ference between the two graphs. It has been shown that 
the latter contribution is zero in the dilute-dense frame¬ 
work |59j . Following this guidance, we will compute the 
two contributions to Eqj^ This issue is further discussed 
in Sec. 4.3.1 of the recent review article [60] . 

The quantity of experimental interest is the two parti¬ 


cle correlation of charged hadrons that is given by 
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where ?7inin,max denote the pseudo-rapidity range of the 
detector and qmm.max (beiiote the transverse 

momentum ranges of trigger and associated particles. 
The function A {ijp, rjq) takes into account the acceptance 
of the uncorrelated background for direct comparison to 
experimental measurements. The exact form of A {rjp, rjq) 
for the GMS, ATLAS and ALICE collaborations has been 
studied in detail and can be found in Ref [TT]. As in 
the previous analyses the fragmentation function D{z) is 
taken to be the NLO KKP parametrization [^ . The 
AT-factor in Eq. j^was taken as AT = 1.5 as in our previ¬ 
ous analysis of p-fp collisions at ^/s = 7 TeV in order to 
accommodate higher order corrections, background un¬ 
certainties and sensitivity to the choice of fragmentation 
functions. 

The number of reconstructed charged hadron tracks 
defined as 





is estimated from the single inclusive gluon distribution 


dA^ 
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J $A(k,j,)4’B(PT -k,r). (7) 


In the above expression is the minimum range of 
momentum for charged particle reconstruction in the ex¬ 
perimental measurement. The ATLAS and CMS collab¬ 
orations use p™“ = 0.4 GeV, which indicates that A(5^“ 
is dominated by soft processes whose uncertainties 
are absorbed into the nonperturbative constant Kg that 
represents a gluon to hadron conversion factor. Because 
Eq. [^ is proportional to 5 '_l , the overall nonperturbative 
prefactor of Eq. [6]is Kg x S± and is fixed to reproduce the 
experimental value of A)I“ for the min-bias Qg = 0.168 
GeV^ multiplicity. 

The variation of A(I“ only comes from the choice of Q§ 
in our calculation, and as noted earlier, high multiplicity 
events can be modeled by using larger values of Qq. We 
find the highest multiplicity bins analysed by CMS and 
ATLAS collaborations correspond to values of Qq that 
are 5-6 times its min-bias value. This is consistent with 
the large values of Qq needed to explain the broad tail 
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FIG. 1. Left: The near side integrated yield Tint versus multiplicity at ^/s = 7 and 13 TeV. Right: Tint versus at ^/s = 7 

and 13 TeV. The experimental data is for the multiplicity bin > 110 mnE]. The lower theory curves have ~ 100 

corresponding to Qq — 0.840 GeV^ at 7 TeV and Qq — 0.672 GeV^ at 13 TeV. The upper theory curves have ~ 125 

corresponding to Qo = 1.008 GeV^ at 7 TeV and Qo = 0.840 GeV^ at 13 TeV. 


of the of the multiplicity distributions in [51] : indeed as 
shown there, large fluctuations of the proton saturation 
scales are essential for a quantitative description of high 
multiplicity events in p+p collisions. In the case of p+Pb 
collisions, the scenario is more complicated since there 
are two saturation scales involved in the problem |64| . 

The strength of the di-hadron correlation defined in 
Eq.[5] is normalized by the number of trigger particles 
defined as 


.^trig — 


dr] / d^Pj. / dz 


Djz) dN /pT\ .. 
z 2 drjd^p^ \ z J ’ ^ ’ 


achieved for a wide range of p-bp 7 TeV data without the 
need for fine tuning. We will now extend our study of 
two particle correlations to the recent ATLAS y/s= 13 
TeV data. The dependence on the center of mass energy 
[65] enters through the parton momentum fraction x ap¬ 
pearing in Eqs. Si and In the CGC framework, it 
is that determines both single and double inclusive 
production in p-|-p collisions. Thus in both and 
the dependence on the center of mass energy enters only 
through Q^, 

= ao) 


In this study, we restrict our comparison with experi¬ 
mental data to p,j, > 1 GeV. For p^ < we expect 
corrections from multiple scattering to be important and 
therefore do not show any results below I GeV. The 
two particle correlation has been computed only recently 
for all pT’s using numerical classical Yang-Mills simula¬ 
tions [5B]. However, these results, which are obtained 
for gluon correlations, cannot be compared directly to 
data in the absence of a satisfactory implementation of a 
hadronization scheme. 

The free parameter S± in Eq. [^does not appear in the 
quantity of experimental interest-the integrated associ¬ 
ated yield per trigger, defined to be 


Y; 


int 




(9) 


where = Ac/f^in. corresponds to the minimum of the 
di-hadron correlation and the second term in the paren¬ 
thesis of Eq. appears due to the ZYAM (zero yield at 
minimum) subtraction employed by the experiments. 

With the above mentioned framework, and a few 
fully constrained parameters, remarkable agreement was 


where = Qg(Qo,i/s). Fixing the value of 
fixes Qg and therefore the value of I1nt(Q|)- This leads 
to a universal scaling curve when Yjnt is plotted against 
at different energies as demonstrated in the left plot 
of Fig. On the same plot we show comparisons to 
recent preliminary data form the ATLAS collaboration 
at 13 TeV |3I] along with the previously measured GMS 
data at 7 TeV mm- Our framework reproduces very 
well this scaling feature of the data. We hope similar 
measurements in p-bp collisions at RHIC will enable one 
to test this prediction at lower energies albeit, it must 
be noted, the larger mean x values at the lower energies 
may be sensitive to physics outside the regime of validity 
of this framework. 

The experimental data points shown in Fig. for 
^ 5=13 TeV are obtained for an acceptance of 2 < | A??] < 5 
and —Aiprain < A(()< Ai^min whereas the results for y/s=7 
TeV is obtained for 2< |A? 7 | <4 and 0< A(^< A^min [3T] . 
The range of transverse momentum for associated and 
trigger particles is same for both energies and chosen to 
be I < < 2 GeV. In order to take into account 

the acceptance difference in A 77 and A(j) the values of 17nt 
at y/s=7 TeV are scaled by a constant factor of 3.6. We 
follow the same approach for the calculation of l^nt in 
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FIG. 2. The variation of odd (left) and even (right) harmonic coefficients an defined in Eq |12| extracted from the A(j) distribution 
of di-hadron yield Y (A0) from ATLAS p+p data |31| at ^ = 13 TeV over 2< |A? 7 | <5. The window of trigger and associated 
transverse momentum 0.5 < < 5 GeV. The a„s are plotted against the reconstructed charged tracks The point 

corresponding to = 130 is measured for the multiplicity bin > 120. (See text for details.) 


our framework for direct comparison with the data. 

We match the value of the at y/s = 7 TeV in our 
computation for Qq = 0.168 GeV^ to the min-bias value 
of CMS data which is = 14. We use this normaliza¬ 
tion to obtain higher by increasing the value of Q§ 

in integer multiples of 0.168 GeV^ for both energies stud¬ 
ied here. Over the range of shown in Fig. the six 
points in our calculation for ^/s = 7 TeV correspond to 
values of Qq in the range of 0.168-1.008 GeV^ and the five 
points for ^/s = 13 TeV correspond to a variation of Qq 
in the range of 0.168 — 0.840 GeV^. The measurements 
of for both CMS and ATLAS do not include effi¬ 
ciency corrections. Since we use the same normalization 
factor for comparison to both CMS and ATLAS data, 
the relative difference in the efficiencies for these two ex¬ 
periments gives rise to about a 5% uncertainty in the 
normalization of our with the experimental data. 

Furthermore, it is not clear that the min. bias configu¬ 
ration of the proton obtained from HERA DIS data is 
the same as the min. bias conhguration of the proton in 
hadronic collisions. There is thefore a considerable sys¬ 
tematic uncertainty (which we estimate |66) to be on the 
order of 10%) of the scaling of the x-axis in the left plot of 
Fig.E We avoid such fine tuning since it does not affect 
the qualitative features of the results presented here. 

In the right plot of Fig. we show the variation of 
Vint with the transverse momentum of trigger and asso¬ 
ciated particles for a fixed bin of > 110. This bin of 
multiplicity corresponds to 0.672 < Qq < 0.840 GeV^ for 
Vs = 13 TeV and 0.840 < Qg < 1.008 GeV^ for = 
7 TeV in our calculation. The values of Vnt correspond¬ 
ing to these two values of Q§ give rise to the two sets 
of curves shown here. Within the regime of validity of 
our calculation (p)^®'=4rg > j GeV), the agreement with 
the data is very good. The two sets of results shown in 
Fig.[T] indicate that for p-l-p collisions Vnt only depends 
on As noted previously, a natural explanation is 

that, for symmetric systems, multiparticle production is 


determined by a single scale Q^. 

In asymmetric systems like p-l-Pb two different satu¬ 
ration scales control multiparticle production. Since the 
evolution of the proton scale with energy is different from 
that of nuclei, a complicated relation between the two 
saturation scales will determine the correlation of 
and Vnt- It will therefore be very interesting to see the 
same measurement of versus Vnt in p-fPb collisions 
at two different center of mass energies. This may be fea¬ 
sible with forthcoming data from 200 GeV p+A collisions 
at the RHIC and in the future at the highest p-l-A ener¬ 
gies of the LHC. A detailed calculation in our framework 
is left for future work. 

We will consider now another important feature of 
the p+p data at -^/s =13 TeV measured by the AT¬ 
LAS Collaboration. In Ref [3T], the Acj) dependence of 
V(At/)) is provided in different bins of multiplicity for 
0.5 < ^ 5 GeV. The full At/) distribution contains 

a lot more information than the integrated Vint • We per¬ 
formed a decomposition of the azimuthal structure of the 
two particle correlation by fitting the experimental data 
to the functional form 

10 

Y{A(I)) = Go+ '^22anCOsnA(j). ( 12 ) 

n—1 

In Fig. [^we show the results of fitting the ATLAS data 
by plotting the different a^s as a function of The 

Fourier components of the two particle correlation show 
a remarkable pattern; a„ < 0 for the odd harmonics and 
On ^ 0 for the even harmonics up to the highest multi¬ 
plicity window of >120. The oq contribution is the 
uncollimated background per trigger and scales linearly 
with AT//'^-this is natural since this corresponds to the 
uncollimated double inclusive yield scaling as 
The negative value of oi is also to be anticipated since 
it corresponds to the azimuthal correlation of mini-jets. 
The 02 contribution is the ridge and is clearly positive-in 
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a hydrodynamic scenario, its square root would approx¬ 
imately correspond to the root-mean-square elliptic flow 
coefficient V 2 - Most interestingly, the 03 contribution is 
negative. The consequences of this can be understood 
by the relationship between the a„’s and the t>„’s as ob¬ 
tained from the single particle distributions with respect 
to a reaction plane. 


(Pt 7 Upi Qt ’ Vq ) 

^o(Pt j Upi Qt ; Uq) 


= Wn(PT)2/p)^n(qT.2/g) • 


(13) 


From the above, it is very difficult to acquire a nega¬ 
tive 03 in hydrodynamics as it is a product of U 3 ’s at 
different rapidities. For a boost invariant expansion, 03 
is necessarily positive definite. Indeed, this factorization 
is observed in Pb-|-Pb collisions and provides strong ev¬ 
idence for collective behavior in nucleus-nucleus (A-A) 
collisions. The even-odd pattern continues for 04 , 5 ^ 6 “ 
these are clearly very small but appear distinct within 
the impressive statistical accuracy of the ATLAS mea¬ 
surement. 

A two particle correlation measurement in hydro¬ 
dynamics (which requires event-by-event simulations) 
would require all o„’s to be positive definite. The AT¬ 
LAS measurement may therefore appear to conclusively 
rule out a hydrodynamic interpretation of the high mul¬ 
tiplicity p-t-p data. However if one construes 03 as result¬ 
ing from a “mini-jet” mechanism plus a “hydro medium” 
mechanism, it is conceivable that a relatively large nega¬ 
tive 03 from the mini-jets might mask a positive 03 from 
the hydro mechanism which would grow with if the 
mini-jet contribution is nearly constant. But such a dis¬ 
tinction between mini-jet and medium is completely ar- 
tihcial and requires fine tuning. Particles of a given px 
do not carry labels as belonging to the one or the other 
particle production mechanism and rescatterings should 
be expected to affect both equally. 

A consistent treatment of both types of two particle 
correlations is available within the CGC/glasma frame¬ 
work itself, which also provides the mechanism, in large 
systems, for matter to thermalize and flow hydrodynam- 
ically. The glasma graph approximation to the full for¬ 
malism is a good approximation when is not too 

large. In this approximation, the anticipated ordering of 
the ttn moments is precisely along the lines of what one 
sees in Fig. 


However when becomes large, coherent multiple 
scattering effects within domains of size 1 /Qs should be¬ 
come important, requiring numerical solutions of Yang- 
Mills equations that account for such effects. In [55], it 
was demonstrated that they can lead to a positive 03 . It 
has not yet been shown that this 03 dominates the one 
from mini-jets computed in the same framework. The 
appearance of a positive 03 in p-l-A collisions, in contrast 
to p-bp collisions at the same may indicate plausi¬ 
bly that this coherent multiple scattering mechanism is 
precocious in p-bA collisions. 

Despite the many caveats noted here, and earlier in 
the text, we believe it is fair to say the high multiplicity 
p-bp results challenge the hydrodynamic interpretation of 
the p-bp ridge. In contrast, the data are consistent qual¬ 
itatively and even quantitatively with the expectations 
of the glasma initial state framework. Further clarity 
can be achieved once high multiplicity data from RHIC 
becomes available. A finer binning in px of the a„ mo¬ 
ments, the extraction of the two particle yields to higher 
values of and the possible extraction of 4-particle 
cumulants would also be helpful. While these are ex¬ 
perimentally difhcult to extract in p-bp collisions, they 
would add clarity to the ongoing debate on the nature of 
the strongly correlated QCD matter produced in small 
systems. 
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